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Abstract 



We discuss the non-linear growth of the kurtosis of the smoothed peculiar velocity field (along an 
arbitrary direction), in an Einstein-de Sitter universe, induced by Gaussian primordial density 
fluctuations. Applying the perturbative theory, we show that, for different cosmological models, a 
departure from the original Gaussian distribution is gravitationally induced only on small scales 
(< 20 — 30 Mpc). Models with scale-free power spectrum P(k) oc k n , with —1 < n < 1, are 
also considered. When the fluid particles move according to the Zel'dovich approximation the 
probability distribution of the peculiar velocity field remains unaltered during the evolution. 

Subject headings: Cosmology - Galaxies: clustering - large-scale structure of the Universe 



1 Introduction 

In the past decade astronomers and cosmologists devoted a great effort to measure large-scale de- 
viations from the Hubble flow and to interpret the cosmological inferences. Redshift-independent 
distance estimators and all-sky catalogs of galaxy redshifts allowed to estimate the (radial) pe- 
culiar motions of galaxies. Undertaking detailed maps of the peculiar velocity field allows for 
example to make dynamical estimates of the density parameter Q a , to test the primordial density 
fluctuation power spectrum P(k), and to measure directly the underlying total (luminous plus 
dark) mass distribution, once the gravitational instability picture is assumed. 

The statistical properties of the large-scale motions may be described in terms of the peculiar 
velocity correlation tensor (Gorski et al. 1989; Groth, Juszkievicz, & Ostriker 1989; Tormen et 
al. 1993) and/or the central moments of the velocity probability distribution function (pdf) p(v) 
(Kofman et al. 1994; Bernardeau 1994). Assuming that the very early density pdf is Gaussian, it 
results that, during the linear regime, p(v) is a Gaussian too, with mean (v) = and variance c 2 = 
(.£f 2 f2 L2 /27r 2 ) Jq°° dkP{k) , H being the Hubble constant. Even if the primordial density field is 
Gaussian distributed, the non-linear time evolution will ensure that the mass density fluctuations 

5 become highly non-Gaussian (Peebles 1980; Fry 1984; Goroff et al. 1986; Juszkiewicz, Bouchet, 
&: Colombi 1993; Catelan & Moscardini 1994), implying a modification of the original Gaussian 
pdf p(v). However, due to the isotropy of the cosmological velocity field, it is expected that all 
the odd moments of p(v) remain zero during the growth of the density fluctuations. A particular 
case is the third central moment, the velocity skewness, discussed for example in Ruamsuwan h 
Fry (1992). Thus, any gravitationally induced departure from the Gaussian distribution may be 
sought only in the fourth central moment of p(v), the velocity kurtosis, as well as in higher order 
even central moments (Grinstein et al. 1987; Kofman et al. 1994; Bernardeau 1994). These have 
been also analyzed in the framework of the global texture model by Scherrer (1992) and Catelan 

6 Scherrer (1994). 

In this work, we study the non-Gaussian content of the velocity pdf in terms of the velocity 
kurtosis as induced by the gravitational growth of the initially Gaussian density fluctuations. The 
velocity kurtosis describes important features such as sharpness of the velocity pdf and the extent 
of its rare-event tail. Taking advantage of the exact perturbative technique (Fry 1984; Goroff et 
al. 1986) and the Zel'dovich approximation (see Grinstein k, Wise 1987), we estimate the kurtosis 
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of the peculiar velocity along a direction a, namely the parameter K v = [(v£) — 3 (i> 2 ) 2 ]/ '(v 2 ) 2 , 
after smoothing with a Gaussian filter. In particular we study the dependence on both the scale 
in the context of several cosmological scenarios and the primordial spectral index, when an initial 
(scale-free) power spectrum P(k) oc k n is assumed. 

The layout of this paper is the following. In Section 2 the exact perturbative theory and 
the (Eulerian version of the) Zel'dovich approximation are reviewed. In Section 3 we discuss the 
gravitationally induced velocity kurtosis parameter K v of an initially Gaussian peculiar velocity 
field. Our results and conclusions are presented in Section 4. 

2 Non— Linear Time Evolution 

We assume that present-day structures formed by gravitational instability from Gaussian fluc- 
tuations 5 in a pressureless fluid with matter density p = pb[l + 5], where pb is the back- 
ground mean density. The density fluctuation field 5 may be written as a Fourier integral, 
5(x,i) = (27r)~ 3 / dkd(k, t) e lk ' x , where x and k are the comoving Eulerian coordinate and 
wavevector, t is the cosmic time. The power spectrum P(k) fully determines the statistics of 
the primordial Gaussian density field, whose variance is a 2 = (l/2ir 2 ) dk k 2 P(k). We filter 
the field 5 by means of a Gaussian window function Wr(x) = (2-kR 2 )^/ 2 exp(— x 2 /2R 2 ) . The 
mass variance on scale R, a\ , is related to P(k) by ojj = (1/2-7T 2 ) Jq°° dk k 2 P{k) [Wr(/c)] 2 , where 
Wji(k) is the Fourier transform of Wr(x) . 

2.1 Equations of Motion: Perturbative Theory 

The time evolution equations for the matter density fluctuation <5(x, t) and the peculiar velocity 
field v(x, t) are the Euler equation and the continuity equation, i.e. 

<9 v+-(v-V)v + ^v = g, (1) 

d o 5 + -V-(l + 5)v = 0. (2) 

a 

Here d = d/dt and spatial derivatives are with respect to x. The density contrast 5 is related 
to the Newtonian gravitational potential A(x, t) = — (47r) _1 / dx! <5(x', t)/\x! — x| via the Poisson 
equation, V 2 A = S . In terms of A the peculiar gravitational acceleration is defined as g = 
— 4ttG pb a VA . We analyze these equations assuming an Einstein-de Sitter universe with no 
cosmological constant. In such a model, the scale factor a is proportional to i 2 / 3 during the 
matter dominated epoch, and the adiabatic expansion implies that QttG pb t 2 = 1 . 

The first-order solution for 5 has the well-known self-similar form, namely, considering only 
the growing mode, S^\x, t) = D(t) <5i(x) , where D(t) oc a(t) is the time growth factor of the mass 
fluctuations. In the linear regime, the peculiar velocity field is proportional to the gravitational 
acceleration 

v = -a§VA. (3) 

This relation shows that the linear velocity field is irrotational; its growing mode corresponds to 
the growing mode of the density field and, for a flat universe, the classical law v = g t ~ t 1//3 is 
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recovered. According to Eq.(3), the particles of the gravitating fluid move along the direction of 
the gravitational force. 

It may be useful to give explicitly the Fourier transform of Eq.(3), i.e., for the component 
along a fixed direction 2, 

D k a ~ 

v a (k,t) = ia— p S(k,t) . (4) 

Higher order approximations of the density solution may be recovered if one expands the mass 
density fluctuation field 5(x, t) about the background solution 5 = 0, namely 5 = J2 n <^"^ with 
S( n ) = 0{5[ l ), then solving the differential equation for any <5( n ) (Peebles 1980; Fry 1984). The 
perturbative expansion for 5 reads (e.g. Goroff et al. 1986): <5(x, t) = Y^=i[D{t)] n <5„(x) . The 
first term of the expansion corresponds to the linear approximation. We see that the scale factor 
D{t) acts as a coupling constant, since 5^ oc D n . 

In a similar fashion, expanding v about the solution v = 0, one obtains 

v(x,f) = -^[D(rv tl (x), (5) 

n=l 

where v„ = 0(v" ) and — V • vi = 5\ . We assume that at any order V A v ra = (Kelvin 
circulation theorem). 

Here we review the exact perturbative technique to solve approximately the equations of 
motion (1) and (2) explicitly up to third order in the peculiar velocity field. In particular, we 
use the third-order solution to compute the fourth-order moment (namely the kurtosis K v ). We 
adopt the same notation of Fry (1984) and Catelan & Moscardini (1994). 

(i) Second-Order Velocity Solution 
The second-order peculiar velocity v^ 2 ) is a solution of the equations 

9o ( av (2)) + ( v (i). V ) v ( 1 )=ag( 2 ), (6) 
5o5 (2) +a -i V .( v (2) + ^(D v (i)) =0) (?) 

where the second-order peculiar acceleration is g( 2 ) = — AttG pb a VA^ 2 ) and V 2 A^ 2 ) = 5^ 2 \ The 
second-order density contribution 5^ has been derived by Peebles (1980). Since 5^ cx D 2 , it 
results that g( 2 ) oc PbaD 2 and the second-order velocity solution reads 



(2) D 
V =~ a D 



2VA( 2 )-5«VA( 1 ) 



+ F 2 , (8) 



where F 2 is a divergenceless vector such that V A v^ 2 ) = (Catelan et al. 1994). The Fourier 

(2) 

transformed v a may be directly obtained from the solution (8), 

^ ) (Kt)=^a^ J ^ i^ 2 )(k',k-k') ^)(k',tH«(k-k',t) , (9) 
where we have defined the kernel 
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It is straightforward to show that (v( 2 )) =0. In an Einstein-de Sitter universe v^ 2 ) ~ t, 
and it grows slower than 5^ ~ i 4 / 3 . We stress the fact that v< 2 > is not parallel to the second- 
order acceleration [oc — VA^ 2 ) ]: this is a consequence of non-locality. Thus, the gravitational field 
changes direction and the particles are not accelerated in a fixed direction, unlike the linear regime. 
The density-velocity relation in the quasi-linear regime and the cosmological implications of non- 
locality have been recently explored by Nusser et al. (1991) and Gramann (1993a). The solution 
(9) is also derived by Gramann (1993b), who applies a second-order Lagrangian perturbative 
technique. 

(ii) Third-Order Velocity Solution 
The third-order approximation is a solution of the differential equations 

d (av< 3 )) + (v« • V) v< 2 ) + (V 2 ) • V) v« = ag( 3 ' , (11) 

do £(3) + a -i V . (V 3 ) + <5« v( 2 ) + 8^ v«) = . (12) 

Here g( 3 ) = — AttG pb a VA^ 3 ) and V 2 A^ 3 ) = 5^ 3 \ The third-order density solution has been 
obtained by Pry (1984). Since g( 3 ) oc pt>aD 3 , and using the results of the previous subsection, 
the velocity v^ 3 ) may be written as 



»« = -.§ 



3 VA( 3 ) + <5« 2 VA« - 2 SV> VA( 2 ) - 5^ V5« 1 + F 3 . (13) 



Again the additive term F3 is such that V A = 0. The Fourier transform of the previous 
expression is 

v®(k,t) = 

J dki y 3 fe(Ek ft -k) ^ 3 )(k 1 ,k 2 ,k 3 )^)(k 1 ,t)^)(k 2 ,t) ( 5( 1 )(k3,t), (14) 

^ ' h=l 

where the third-order kernel is 

K( 3 )( kl ,k 2 ,k 3 ) = 3 J( 3 )( kl ,k 2 ,k 3 ) - ^1 J( 2 )(k 2 ,k 3 ) - ^l^y ^ (2) (kl,k 2 ) , (15) 

and the functions and j( 3 ) , corresponding to the second- and third-order density solutions 
(see e.g. Fry 1984; Catelan & Moscardini 1994), read respectively 

7(2) /T, I, ^ 5 , ki-k 2 , 2 ^ kx-k 2 ^ 2 



C 2 I \ K\ K2 

j( 3 )(kx,k 2 ,k 3 ) = J( 2 )(k 2 ,k 3 



1 1 ki • (k 2 + k 3 ) 4 k • ki k • (k 2 + k 3 ) 



3 3 (k 2 + k 3 ) 2 9 kl (k 2 + k 3 ) 2 

2 k • k, k (k 2 + k 3 ) (k 2 + k 3 ) • k 3 , 1 k • k 2 k ■ k 3 

.2 
3 



9 kl (k 2 + k 3 ) 2 kl + 9 ^ ki • (17) 



It is not difficult to show that (v^ 3 )) = 0. In an Einstein-de Sitter universe, ~ 
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(Hi) General Solution 

As Goroff et al. (1986) have shown, the (Fourier transformed) n-th order velocity solution may 
be represented in integral form as 

^( k ) = *p|ll J J^~ d ^j [(2^D(E k ^- k )]^ (n) (ki,-..,k n ). (18) 

The presence of the Dirac delta function comes from momentum conservation in Fourier space. 
The kernels are homogeneous (with degree 0) functions of the wavevectors ki, . . . ,k n , de- 
scribing the effects of non-linear collapse (tidal and shear effects). In general the are very 
complicated for n > 3. [A discussion of the properties of the kernels is given in Wise (1988). 
Explicit recursion relations with their Feynman diagrammatic representations are given by Goroff 
et al. (1986) and Wise (1988).] 

2.2 Zel'dovich Approximation 

In the Zel'dovich approximation (Zel'dovich 1970) the motion of particles from the initial comoving 
(Lagrangian) positions q is approximated by straight paths. The Eulerian position at time t is 
then given by the uniform motion 

x(q,t)=q + D(*)S(q), (19) 

where D(t) is the growth factor of linear density perturbations and S(q) is the displacement vector 
related to the primordial velocity field. Grinstein & Wise (1987) give an Eulerian representation 
of the Zel'dovich approximation by a diagrammatic perturbative approach similar to that of the 
previous section. They showed that the n-th order perturbative corrections <5 n (x), when the 
density fluctuation field 5 is evolved according to the Zel'dovich approximation, are such that 

«(x, *) = E4r- W)]» s fcl ... ShB . (20) 

n=l n - { hn ]=l° Xh i ° Xh " 

Here J2[h„] = ^2hi ' ' '^2h n - Note that the first term recovers the linear approximation, in that 
S = vi, where <5i(x) = — V • Vi . This expansion for 5 corresponds to different symmetric kernels 
K-ZA ' wn i cn can be written in the following compact form 



4i(ki,...,k n) =in^, (2D 

n - h=l K h 



where k = J2h=i kfe • The kernels are the same obtained in the expansion of the density 

contrast (see Catelan & Moscardini 1994) and are symmetric by construction. 

3 Kurtosis of the Velocity Field 

In this section, we compute the gravitationally induced kurtosis K v of an initial Gaussian velocity 
field in a flat universe. We restrict the calculation to the velocity along a chosen direction a . The 
lowest order non-zero reduced contribution to K v is 

^ 1)2 ) 2 ^-6(,« 2 ,( 2 ) 2 ) c + 4(^) 3 ,( 3 )) c , (22) 
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where the subscript c indicates the connected part of the four-point velocity correlation. 

It is not difficult to verify that K v depends on the normalization of the power spectrum. 
Furthermore it varies with time. The related quantity which is independent both of the power 
spectrum normalization and of time is S^ v = K v / (8^ 2 ): due to these properties, observational 
estimates of S^ v would allow one to detect intrinsic features of the evolved large-scale velocity 
field. 

From the perturbative results in Eqs.(9), (14) and (22), one finally gets the integral expression 
of the kurtosis of the smoothed velocity field vr: 

• / m __24_ / b\ 4 f rfki dk 2 dk 3 kf k% Aff (- kf - fcg - Aff) 
V{) ~<^\D) J ~ J?nf fcf fc|fc||k 1+ k 2 + k3| 2 X 

W R (h) W R {k 2 ) W R (h) W^(|ki + k 2 + k 3 |) x 

P(k 3 ) Af )(k l5 k 2 , k 3 ) + 2 K^(-k 2 , k 2 + k 3 ) Af )(k 1; k 2 + k 3 ) P(|k 2 + k 3 |)] , 

(23) 

where a 2 v R is the variance of smoothed velocity field. The kernels are obtained by complete 
symmetrization of the kernels K^ n \ The analogous expression in the Zel'dovich approximation 
is obtained by replacing the kernels with the corresponding . Finally, note that P(k) 
completely describes the process of growth of the higher order velocity moments from Gaussian 
initial conditions. 



P(k 1 )P(k 2 ) 



4 Discussion and Conclusions 

We calculate the previous integrals, by an Adaptive Multidimensional Monte Carlo Integration 
subroutine, in the framework of several different cosmological models. In particular, we consider: 

(1) standard cold dark matter model (SCDM), i.e. with 6=1 for the linear biasing parameter; 

(2) a biased (BCDM) version of the same model, with b = 1.5; (3) tilted cold dark matter model 
(TCDM), with spectral index n = 0.7 and b = 2.0; (4) a 'mixed' model (MDM), with 60% cold 
and 30% hot dark matter, b = 1.5. All transfer functions have been taken from Holtzmann (1989), 
with the Hubble constant H Q = 50 km s _1 Mpc -1 . 

The results for K v are shown in Fig. la. The associate relative uncertainty estimates from the 
Monte Carlo Integration (not shown in the figure for clarity) are always smaller than 5%. For 
any model the velocity kurtosis is a decreasing function of the smoothing scale R. On large scales 
K v is always consistent with zero, thus the distribution of the peculiar velocity field is very close 
to a Gaussian one. However, a strong non-Gaussian signature - i.e. K v > 1 - is induced on 
scales smaller than < 20 — 30 Mpc, where nevertheless observations are at the moment affected 
by large uncertainties. Furthermore, even if the general trend is similar for all considered models, 
the departure from Gaussianity is larger for models with low biasing parameter: TCDM gives 
always the lowest values of K v . 

A previous estimate of the kurtosis of the velocity field was obtained by Kofman et al. (1994), 
who used the smoothed velocity fields of a A^-body simulation of the standard cold dark matter 
model with box-size of 400 Mpc to study the velocity distribution p(v). In order to have a more 
direct comparison, in this case we calculated K v using their choice for the CDM transfer function 
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(Davis et al. 1985) and limiting the numerical integration of the expression (23) in the range of 
wavevectors spanned by the iV-body simulation: from the frequency corresponding to the box- 
size to the Nyquist one. In Fig. lb we show both our perturbative evaluations and the Kofman et 
al.'s data: the agreement is quite good, except for the small scale (12 Mpc) result. It might seem 
surprising that the numerical simulations, which are supposed to better describe the highly nonlin- 
ear regime, miss the departure from Gaussianity of the velocity field: in fact, for the density field, 
the perturbative evaluations of higher moments underestimate the iV-body results. A possible 
explanation is that the multiple filtering (a trilinear interpolation plus a small-scale smoothing) 
necessary in the simulations to reconstruct the velocity field from the particle distribution might 
smooth the small-scale non-Gaussian signal. Moreover, a further "numerical smoothing" can 
appear due to the intrinsic gridding present in the particle-mesh code. In any case, it is necessary 
to be cautious in the measurements of higher order moments in iV-body simulations because they 
are strongly affected by the high-value tails, which characterize the particular realization, size of 
the box and other numerical problems (Kofman, private communication): this is also shown by 
the large error bars in the iV-body results. 

We also calculate the velocity kurtosis in terms of the intrinsic parameter S^ v . In Fig. 2a, £4 „ 
is plotted for the same models previously considered: in this case, due to the independence of 
v on the power spectrum normalization, the biased and the standard CDM originate the same 
curve. All models present a similar behavior and the differences are always inside the error bars, 
shown for clarity only for SCDM. In Fig.2b, the dependence of on the primordial spectral 
index n is shown for scale-free power spectra P(k) oc k n , with n in the range —1 < n < 1, for 
both the perturbative and Zel'dovich approximations. Due to the assumed scale-invar iance, S^ v 
only depends on the primordial spectral index n, and not on the scale R. It may be noted that 
in the Zel'dovich approximation, unlike the perturbative case, the velocity kurtosis is practically 
constant and consistent with zero for any value of n. Thus, we confirm the results of Kofman et 
al. (1994), who demonstrate that the Eulerian Gaussian one-point pdf p(v) is time-invariant as 
long as the Zel'dovich approximation holds. This is essentially due to the simple time scaling of 
the particle Eulerian position x(q, t) in Eq.(19). 

Finally we want to stress that one has to be careful about making quantitative comparisons 
between our results, directly related to the underlying (dark plus luminous) mass distribution, and 
observational data. In particular, Grinstein et al. (1987) suggest that when point-like luminous 
objects, like galaxies, are used to sample the peculiar flow within a region of the sky, it is the 
volume average of n(x)v(x) which is actually measured instead of the volume average of v(x), n(x) 
being the number density of luminous tracers. If the objects are biased tracers of the underlying 
mass distribution, then nonlinear effects on small-scales may preclude any direct comparison of 
the observed velocity moments with ensemble expectations. A further analysis of this problem is 
in progress. 
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Figure captions 



Figure 1. The kurtosis ratio K v of the velocity field. Left panel: the behavior versus the scale R 
for Gaussian filter and different cosmological models (parameters are in the text): standard cold 
dark matter (solid line); biased cold dark matter (dotted line); tilted cold dark matter (dashed 
line); mixed dark matter (dotted-dashed line). Right panel: comparison with results from N— 
body (Kofman et al. 1994, filled squares) in the case of standard cold dark matter. Filled triangles 
and solid line refer to the perturbative estimates when the integration is limited in the interval 
ranging from the box-size frequency to the Nyquist one; dotted line refers to the perturbative 
estimates when the integration is over all frequencies. 

Figure 2. The intrinsic parameter Si tV . Left panel: the behavior versus the scale R for different 
cosmological models: standard cold dark matter (triangles), tilted cold dark matter (squares), 
mixed dark matter (circles). Error bars, shown only for SCDM, refer to the associated uncertainty 
estimate from the Monte Carlo Integration. Right panel: the behavior versus the primordial 
spectral index n for power-law spectra P(k) cx k n and Gaussian filter, for both the perturbative 
(triangles and solid line) and Zel'dovich approximations (squares and dotted line). 
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